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ABSTRACT 

We examine the process of element segregation by gravity in giant elliptical galax- 
ies and X-ray clusters. We solve the full set of flow equations developed by Burg- 
ers for a multi-component fluid, under the assumption that magnetic fields slow the 
diffusion by a constant factor Fb- Compared to the previous calculations that ne- 
glected the residual heat flow terms, we find that diffusion is faster by ~ 20%. In 
clusters we find that the diffusion changes the local abundance typically by factors of 
I + 0.3(T/10 8 K) x - 5 /F B and 1 + 0.15(T/10 8 K) 15 /F B for helium and heavy elements, 
respectively, where T is the gas temperature. In elliptical galaxies, the correspond- 
ing factors are I + 0.2(T/10 7 K) 15 /F B and I + 0.1(T/10 7 K) l */F B , respectively. If 
the suppression factor Fg is modest, diffusion could significantly affect observational 
properties of hot X-ray clusters and cD galaxies. In particular, diffusion steepens the 
baryon distribution, increases the total X-ray luminosity, and changes the spectrum 
and evolution of stars that form out of the helium-rich gas. Detection of these diffu- 
sion signatures would allow to gauge the significance of the magnetic fields, which also 
inhibit thermal heat conduction as a mechanism for eliminating cooling flows in the 
same environments. 
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^ 1 1 INTRODUCTION 

x : 

' In a multicomponent plasma, gravity acts so as to separate 
different elements. The thermal equilibrium state for the 
number density of particles of mass mi and temperature T 
is proportional to the Boltzmann factor, e~ mi ^^ r ^ kT , where 
(j)(r) is the gravitational potential. Thus in self-gravitating 
objects formed out of a homogeneous mixture of elements, 
the heavier elements (such as helium, metals) sink towards 
the center and the lighter elements (such as hydrogen) 
rise up. Fabian & Pringle (1977) estimated that in high- 
temperature X-ray clusters significant abundance gradients 
might arise for metals. However, their calculation did not al- 
low for the suppression of diffusion by magnetic fields, which 
could play an important role in clusters. Since a charged par- 
ticle is confined to move along magnetic field lines, diffusion 
can proceed uninhibited only on the coherence length-scale 
of the magnetic field. In clusters the scale of magnetic field 
is likely to be small and the diffusion speed could be sup- 
pressed by several orders of magnitude (Chandran & Cowley 
1998), making its effect negligible. In recent years, however, 
it has been suggested that turbulence may suppress trans- 
port processes by only a modest factor of a few (Narayan 
& Medvedev 2001; Malyshkin 2001). Such models are moti- 



vated by the lack of observational evidence for cooling flows 
in the cores of clusters and elliptical galaxies (Mollendi & 
Pizzolato 2001; Fabian et al. 2001). Cooling flows would be 
suppressed if thermal heat conduction by electrons is not 
strongly inhibited by the magnetic fields. 

In this paper we estimate the impact of diffusion on ob- 
servational properties of galaxies and X-ray clusters for dif- 
ferent magnetic suppression factors. Unlike previous analy- 
ses (Fabian & Pringle 1977, Rephaeli 1978, Gilfanov & Syun- 
yaev 1984, Qin & Wu 2000, Chuzhoy & Nusser 2003), we 
solve the full diffusion equations derived by Burgers (1969). 
We find that if the suppression factor is modest, then dif- 
fusion would lead to significant abundance gradients for all 
elements. In extreme cases, the abundance of helium (which 
except for deuterium is the fastest diffusing element) can 
grow above that of hydrogen. Because helium sedimentation 
changes the mean molecular weight, the gas density profile 
steepens and the total luminosity of the cluster increases. 
Since the diffusion speed depends strongly on temperature, 
this effect may contribute to the luminosity-temperature de- 
pendence observed in clusters (Chuzhoy & Nusser 2003). 
The change in the mean molecular weight would also affect 
mass estimates of hot clusters and elliptical galaxies under 
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the assumption of hydrostatic equilibrium (Qin & Wu 2000, 
Chuzhoy & Nusser 2003). In addition to its effect on the 
X-ray emitting gas, the increase in the helium abundance 
could affect the spectrum and evolution of stars that form 
out of this gas. 

Because the diffusion speed in a plasma scales with tem- 
perature as T 3 / 2 , all previous work was restricted to X-ray 
clusters which are the hottest virialized objects. However, we 
find that because of the lower gas mass fraction, diffusion 
would proceed faster in hot cD galaxies with T~2x 10 7 K 
than in clusters with T < 10 8 K. 

The paper is organized as follows. In §2.1 we write the 
general diffusion equations. In §2.2 we evaluate the diffu- 
sion timescales and describe the diffusion in the perturbative 
regime, where the fractional deviations of the abundances 
from uniformity are small. In §3 we describe diffusion in the 
non-linear regime, assuming the NFW mass profile for the 
dark matter (Navarro, Frenk & White 1997). We summarize 
our results in §4. 



2 DIFFUSION EQUATIONS 
2.1 General equations 

Each species of particles s is described by a distribution func- 
tion F s (x,v,t) normalized to unit integral, a mean number 
density n s , an ionic charge q s = Z s e and a particle mass m s . 
The cross-section for Coulomb scattering between particles 
of species s and of species t is given by 



<T rt = 2^e 4 Z s 2 Z t 2 (fc s T)- 2 lnA st , 



(1) 



where lnA st is the Coulomb logarithm. For simplicity, we 
approximate the Coulomb logarithm for all species by a con- 
stant, hiA st = 40, which is characteristic for clusters. The 
friction coefficient between species s and t is 



K Bt = (2/3)n s n t ir s t 



f 2kBTm s m t 



\ m, 

The mean fluid velocity is 
u = (£n s m s u s )/(£n s m s ), 



m s m t \ 
+ m t J 



1/2 



(2) 



(3) 



call & Loeb (1994), are 
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EK s t[(w t - w B ) + 0.6(x st r s - y a tn)], (7) 
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-j/ st [1.6a; s t(r s + r t ) + Y st r s - 4.3x st r t }}, (8) 

where g is the gravitational acceleration, E is the elec- 
tric field, x st = m t /(m s + m t ), y st = m s /(m s + mi) and 
Yet = 3j/ st + 1.3x s tTit/m s . From equations (7) and (8) we 
can obtain a solution for the diffusion velocity and the heat 
flow of each species 
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(10) 



where the coefficients If, I T , J? and J T depend only on the 
composition of the gas. In the presence of magnetic fields, we 
introduce a constant suppression factor, Fg 1 , for all these 
coefficients. The actual suppression factor in reality should 
depend on the geometry and fluctuations of the magnetic 
field which may vary with position; but the lack of detailed 
data on these properties calls for a simple-minded approach 
of the type we adopt. 

If the diffusion is "switched-on" when the gas is in hy- 
drostatic equilibrium, then initially there would be no mass 
flow (u = 0). For each sinking helium ion there must be 
approximately four protons and two electrons that rise up. 
This means that there is a net outflow of particles, which 
leads to a decrease in the total pressure (p — T,n B kBT). For 
a fixed mass density, the drop in gas pressure means that the 
assumption of hydrostatic equilibrium (Vp/p = g) can not 
hold once diffusion starts to operate. In practice the sound 
speed of the gas is much larger than the diffusion speed, i.e. 
the dynamical timescale of the system is much shorter than 
the diffusion timescale. Thus, the gas will always be close 
to hydrostatic equilibrium and the above equations for the 
relative velocities of all species remain valid. To find a small 
mean fluid velocity, we solve the equation of motion, 



where u s is the mean fluid velocity of each species. The 
diffusion velocity of species s is defined as 



du Vp 
dt p 



(11) 



W s = U s — U. 



and its "residual heat flow vector" is (Burgers 1969) 

5 



2k B T 



J F s (v-u)|v-u| 2 dv- 



(4) 



(5) 



If the gas is in hydrostatic equilibrium (i.e. Vp/ p = g), the 
Burgers equations for mass, momentum and energy conser- 
vation in spherical symmetry, as formulated by Thoul, Bah- 



where g = g(r) is the gravitational field dictated by the total 
mass distribution. 

If we wish to evaluate the diffusion velocity of hydrogen 
(s = p) and helium (s = a), which are the dominant compo- 
nents of the cosmic plasma, we may neglect the contribution 
of all other species. Thus, the drift velocity of helium ions is 
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(12) 
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where J p , I a and It are functions of hydrogen mass fraction 
X. The velocity of the protons can be obtained simply from 
the relation w p = w a (X — l)/X. 

The diffusion velocity of heavy elements is determined 
primarily by the drag forces from hydrogen and helium 
(Chuzhoy & Nusser 2003). Helium provides an inward ra- 
dial friction force for the metals, while protons push them 
radially outwards. The resulting velocity can be determined 
by balancing these friction forces. Together with the pres- 
sure gradient force, the gravitational and electric forces for 
the metals increase their diffusion velocity relative to hydro- 
gen by a factor ~ (1 + 4/Z) (Chuzhoy & Nusser 2003). For 
heavy elements such as iron (Z — 26) this correction can 
be ignored. Using equation (7) we obtain the approximation 
for the diffusion velocity of all heavy elements 



(1 - X)w a - 1.2(1 - X)r a - 0.6Xr p 



(13) 



The total heat flow and the temperature evolution are 
given by 



Q 

9(|n 3 fcT) 

dt 



X{n s r s kT), 

S 

1 d(r 2 Q) 
r 2 dr 



(14) 



(15) 



Somewhat surprisingly it turns out that Q and dT /dt are 
generally nonzero even in the absence of a temperature gra- 
dient. As implied by equation (10), the density gradients can 
alone be responsible for the heat transport and eventually 
cause temperature gradients. However, because the coeffi- 
cients J B are several orders of magnitude smaller than J T , 
the effect of the heat flow on the temperature profile is neg- 
ligible in absence of an initial temperature gradient. We also 
find from Burgers equations that the density gradients can 
not produce the heat transport and temperature gradients 
after all elements reach equilibrium distribution. 



2.2 Diffusion in the perturbative regime 

The helium to hydrogen ratio starts from the uniform distri- 
bution set by the big bang. In this section, we derive simple 
analytic results for small deviations from this uniform initial 
state. In the interiors of X-ray clusters and elliptical galax- 
ies, the density changes by a much larger factor than the 
temperature and so | | > | | . Since I B and I T are of 
similar magnitude, equation (12) simplifies to 



T 5/2 <91nn p 



(16) 



Initially the hydrogen fraction X ~ 
/ = Ip + IS = 2.6 x W 7 /F B cm-V 



3/4 everywhere and 
^K" 5 / 2 , where F B is 
the suppression factor due to the magnetic field. Using the 
assumption of approximate hydrostatic equilibrium, we can 
rewrite equation (16) as 



w a ,o 



= 60 km s 



(17) 



10- 



) (io 8 k) Go- 3 cm- 3 ) Fb 



The velocity obtained above is ~ 20% higher than the esti- 
mate by Chuzhoy & Nusser (2003) that neglected the heat 
flow terms and ~ 30% lower than the estimate by Qin & Wu 
(2000) that also neglected electric fields. 



Substituting equation (17) into the continuity equation 



gives 



dn a _ n a f fg \ 1 ( 
dt ~ 12 Gyr Vo.l/ V 



T 



10 8 K 



3/2 



(18) 



where f g is the local gas fraction. Equation (18) implies that 
at the initial stages of diffusion, the local helium density 



grows as e 



t/T 



(1 + t/r), where 



= 12G ^ (rn) (ioSk)" 372 ^- 



(19) 



Similarly the density of hydrogen declines as e */ 3t for 
X = 3/4, so that the relative abundance of helium grows 



4t/3r 



The initial diffusion velocity of metals is 



of 



helium velocity, and so the metallicity grows as ~ e 



2t/3r 



3 DIFFUSION IN GALAXIES AND CLUSTERS 

For the radial distribution of the mass density of the dark 
matter, we assume the NFW profile (Navarro, Frenk & 
White 1997) 



Pd{r) 



Ps 



(cr/r vir )(l + cr/r vir ) 2 



(20) 



where p s — const, c is the concentration parameter and 
r V i r is the virial radius. We make separate calculations for 
c = 4 and c = 15, which are characteristic values for clus- 
ters and elliptical galaxies, respectively (Navarro, Frenk & 
White 1997; Wechsler et al. 2002). We assume that the gas 
is initially isothermal and in hydrostatic equilibrium. Ne- 
glecting the contribution of the baryons to the gravitational 
potential, one obtains the following radial mass profile for 
the gas (Makino, Sasaki & Suto 1998) 



p g (r)=p e v (l + cr/r vir ) 



r l r vir/ cr 



(21) 



where r\ = AttG fim p p B r 2 i r / 'c 2 kT '. Assuming that the gas is at 
the virial temperature, n = 10 and n = 16.5 correspond to 
c = 4 and c = 15 respectively. For clusters, we chose the ratio 
po/Ps so that the gas mass fraction within the virial radius 
is equal to the cosmic value of (fi(,/f2 m ) = 0.17 inferred by 
WMAP (Bennett et al. 2003). In elliptical galaxies, typically 
~ 90% of the baryonic mass is in stars (Brigenti & Mathews 
1997), and so we take the gas fraction ten times lower. We 
assume that the initial abundances are uniform and that the 
temperature remains constant in time as well as in space. 

We solved equations (6)-(8) and (11) for the above ini- 
tial conditions. Because the gas was assumed to be isother- 
mal, the evolution of metallicity and helium abundance (fig. 
1-4) depends mainly on the local gas mass fraction (fig. 5 
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iV 10- 1 10° 

Figure 1. Helium to hydrogen mass density ratio vs radius in 
clusters. The numbers on the curves correspond to the elapsed 
time in Gyr divided by F B X (T/10 8 K)- 3 / 2 . 
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Figure 2. Same as Fig. 1 for elliptical galaxies. 



and 6). In clusters, the gas mass fraction declines monotoni- 
cally with radius initially, accounting for the rise in metallic- 
ity and helium abundance towards the center. In ellipticals, 
the gas fraction profile is more complex, which in turn re- 
sults in non-monotonic abundance profiles. 




Figure 3. Metallicity vs radius in clusters. The plot is normalized 
by the initial metallicity Zq. The numbers on the curves corre- 
spond to the elapsed time in Gyr divided by F B X (T/10 S A") _3 / 2 . 




Figure 5. Gas mass fraction vs radius in clusters. The numbers 
on the curves correspond to the elapsed time in Gyr divided by 
F B X (T/10 8 ^)- 3 / 2 . 

4 DISCUSSION 

Our calculations indicate that in the absence of magnetic 
suppression, the diffusion timescale in hot clusters and cD 
galaxies is comparable to the age of these systems. For a sup- 
pression factor Fb i; 5, the diffusion can be well described in 
the perturbative regime (§2.2) for most clusters and galax- 
ies, except in the central regions of very hot clusters. The 
abundance change depends on the gas mass fraction, and 
characteristic numbers were derived by adopting f g — 0.17 



0.05 




1V 2 10- 1 10° 
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Figure 6. Same as Fig. 5 for elliptical galaxies. 
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and f g = 0.017 for clusters and cD galaxies, respectively. In- 
terior to the virial radius of clusters with an age t ~ 5 Gyr, 
the helium and metal abundance would increase by a factor 
of 1 + 0.3(T/10 8 K) 15 /F B and 1 + 0.15(T/10 8 K) 15 /F B , 
respectively. For elliptical galaxies with an age t ~ 10 Gyr, 
the corresponding factors are 1 + 0.2(T/10 7 K) 15 /F s and 
1 + 0.1(T/10 7 KY^'/Fb, respectively. Although the effect of 
diffusion on the distribution of metals is rather modest, its 
effect on the distribution of helium is more substantial. As 
indicated by Figures 1 and 2, some regions inside the hottest 
X-ray clusters or cD galaxies, show an increase of ~ 50% of 
helium abundance for Fb = 3. 

Detection of the diffusion signatures can be used to 
gauge the effect of magnetic fields which also limit ther- 
mal heat conduction from suppressing cooling flows in the 
same environments. Although the increased helium abun- 
dance can not be observed directly, it leads to a number of 
indirect observable signatures. 

As shown by Fig. 5, diffusion results in a substantial 
steepening of the gas distribution. This would in turn lead 
to an increase in the total X-ray luminosity. Since diffusion 
is fastest in high temperature clusters, this process may con- 
tribute to the discrepancy between the observed luminosity- 
temperature relation L oc T 3 (Mushotzky 1984; Edge & 
Stewart 1991; David et al. 1993) and the form L oc T 2 ex- 
pected from self-similar arguments (Kaiser 1986). 

The mass profile of clusters and elliptical galaxies is usu- 
ally inferred from the observed distribution of hot X-ray gas 
under the assumption of hydrostatic equilibrium. Increasing 
the helium fraction by ~ 50% is equivalent to increasing 
the mean molecular weight by ~ 10%, which in turn results 
in overestimating the mass also by ~ 10%. Thus, in cases 
where alternative methods for mass determination (such as 
gravitational lensing or the Sunyaev-Zel'dovich effect) reach 
better precision, these methods may be used to constrain the 
diffusion velocity and the magnetic suppression coefficient. 

The change in the helium fraction would also affect the 
evolution of stars formed from the enriched gas. D'Antona 
et al. (2002), who analyzed globular clusters, found that 
increasing the initial helium abundance by 20% results in 
the formation of a very blue horizontal branch. It would 
be particularly interesting to search observationally for the 
distinct spectral and evolutionary characteristics of helium- 
rich stars within X-ray clusters and cD galaxies. 
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